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Even  though  the  steady  or  quasi-steady  combust' m  of  isolated  droplets  or  droplet 
sprays  has  been  studied  extensively  and  is  reasonably  well  understood,  the  problem  of 
transient  gas  dynamics  coupled  to  droplet  combustion  is  so  poorly  understood  at  the 
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Aside  from  the  detonation  of  spray  droplet  mixtures,  in  which  droplet  breakup  and 
combustion  are  occurring  under  extremely  transient  conditions,  a  full  descriptive 
theory  of  this  phenomena  will  be  invaluable  in  furthering  our  understanding  of  such 
varied  phenomena  as  liquid  rocket  engine  instabilities,  engine  thrust  transients, 
combustion  in  liquid  propellant  gun  systems,  high  pressure  fuel-air  gun  systems,  spray 
ignition  mechanisms,  and  hybrid  rocket  dynamics. 

This  report  presents  a  theoretical  analysis  of  one  important  phase  of  heterogeneous 
flow,  namely  a  model  for  the  aerodynamic  stripping  of  fuel  from  the  assumed  spherical 
droplets. 
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1 . 0  ABSTRACT 


Even  though  the  steady  or  quasi-steady  combustion  of  isolated 
droplets  or  droplet  sprays  has  been  studied  extensively  and  is 
reasonably  well  understood!  the  problem  of  transient  gas  dynamics 
coupled  to  droplet  combustion  is  so  poorly  understood  at  the 
present  time  that  adequate  models  to  describe  the  phenomena  have 
yet  to  be  proposed.  The  problem  is,  however,  extremely  important 
to  the  technology  of  combustion  devices.  Aside  from  the  detonation 
of  spray  droplet  mixture', ,  in  which  droplet  breakup  and  combustion 
are  occurring  under  extremely  transient  conditions,  a  full  des¬ 
criptive  theory  of  thi3  phenomena  will  be  invaluable  in  further¬ 
ing  our  understanding  of  such  varied  phenomena  as  liquid  rocket 
engine  instabilities,  engine  thrust  transients,  combustion  in 
liquid  propellant  gun  systems,  high  pressure  fuel-air  gun  systems, 
spray  ignition  mechanisms,  and  hybrid  rocket  dynamics. 

This  report  presents  a  theoretical  analysis  of  one  important 
phase  of  heterogeneous  flow,  namely  a  model  for  the  aerodynamic 
stripping  of  fuel  from  the  assumed  spherical  droplets. 


2.0  INTRODUCTION 


One  of  the  most  difficult  problems  encountered  in  the  develop¬ 
ment  of  high-speed  compression-ignition  chambers  has  been  the 
proper  atomization  and  distribution  of  the  fuel  in  the  combustion 
chamber  during  the  extremely  short  time  available.  Rapid  combustion 
of  the  fuel  does  not  take  place  as  soon  as  it  enters  the  combustion 
chamber,  but  a  certain  time,  known  as  the  ignition  lag,  elapses 
during  which  the  temperature  of  the  fuel  is  raised  to  its  auto¬ 
ignition  point  by  the  absorption  of  heat  fiom  the  compressed  air. 

The  rate  of  heat  absorption  by  a  fuel  drop  is  directly  proportional 
to  its  surface  area;  the  rate  of  its  temperature  rise  is  inversely 
proportional  to  its  volume.  Because  the  surface  area  varies  as 
the  square  of  the  diameter,  whereas  the  volume  varies  as  its 
cube,  a  small  drop  will  hav.'  a  shorter  ignition  lag  than  a  larger 
one.  The  time  required  for  the  complete  combustion  of  small 
drops  is  also  less  than  that  for  the  large  ones;  therefore,  the 
smallest  drops  are  the  most  desirable  if  they  can  be  obtained 
without  sacrificing  good  distribution. 

This  report  discusses  a  description  of  one  of  the  physical 
phenomena  for  reducing  the  size  of  the  droplets,  namely  that  of 
aerodynamic  stripping. 
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3.0  DESCRIPTION  OF  THE  PROBLEM 
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burning  rates  are  obtainable  which  are  higher  than  those  possible 
under  the  conditions  of  low  velocity  and  forced  convection  where 
no  disintegration  occurs.  In  order  to  increase  the  rate  of  a 
complete  combustion,  the  stripping  effects  of  the  breaking  up 
process  of  a  fuel  droplet  is  considered  in  this  study.  This 
theoretical  study  on  stripping  effects  of  a  droplet  in  high 
speed  air  is  made  in  order  to  investigate  how  to  properly  model 
such  a  physical  phenomena.  Therefore,  the  formulation  of  the 
problem,  i.e.,  the  governing  equations  and  the  proper  boundary 
conditions  are  given  considerable  detail. 


Basically,  this  is  a  two-phase  viscous  flow  problem. 

Since  the  droplet  is  in  a  high  speed  air  stream,  the  problem  is 
further  simplified  by  application  of  standard  boundary  layer 
theory.  This  means  that  a  high  Reynolds  number  flow  problem 
is  reduced  to  a  two-phase  boundary  layer  flow  problem.  In 
order  to  illustrate  this  problem,  a  sketch  of  the  flow  field  is 
presented  in  the  following  sketch: 
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The  two  phase  boundary  layer  flow  problem  has  been  developed 
for  many  engineering  applications,  such  as  ablation,  spallation, 
aerodynamic  shattering,  etc.  Ranger  and  Nicholls^  studied  the 
aerodynamic  shattering  of  a  liquid  droplet.  Their  work  gives 
experimental  results  and  has  attempted  a  theoretical  description. 
The  theoretical  analysis  in  their  study  is  to  integrate  the 
boundary  layer  equation  by  Karman's  momentum  integral  method. 

Then  an  approximate  solution  of  the  velocity  distribution  of  the 
gas  and  the  liquid  is  posed.  This  is  exactly  the  same  as  G.  I. 
Taylor's  analysis*'2-'  in  a  somewhat  different  problem.  With  the 
assumption  of  knowing  a-priori  the  velocity  distribution,  they 
calculated  the  mass  stripped  away  from  the  droplet.  G.  W.  Sutton v 
studied  a  related  problem  of  hydrodynamics  and  heat  conduction 
of  a  melting  surface  near  the  stagnation  point  of  a  high  speed 
stream.  He  carefully  studied  the  boundary  conditions  and  for¬ 
mulated  the  problem  as  an  axially  symmetric  Falkner-Skan  flow. 

He  also  applied  his  analysis  to  experiments  on  the  ablation  of 
reinforced  plastics  in  supersonic  flowv  \  E.  W.  Adams  inves¬ 
tigated  transient  and  quasi-steady  performance  of  melting  type 
re-entry  shield*  .  There  are  many  others  who  also  contributed 
in  this  field;  the  references  are  included. 
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4.0  ASSUMPTIONS  REQUIRED 


Because  this  problem  is  quite  complicated,  some  assumptions 
have  to  be  made  in  order  to  simplify  the  problem,  so  that  it  can 
be  made  mathematically  tractable.  The  assumptions  made  are  as 
follows: 

1.  The  flow  is  considered  as  i  compressible. 

2.  It  is  assumed  that  the  problem  is  one  of  a  steady  state 
flow. 

3.  There  is  no  slip  at  the  interface  of  gas  and  liquid  boundary 
layers  (  at  y  *  0). 

4.  The  outer  solution  of  velocity  over  the  spherical  droplet 
is  assumed  to  be  "Uoq*  ^  sim  ,  which  is  the  poten¬ 
tial  flow,  where  R  is  the  radius  of  the  droplet,  and  x 

is  the  curvilinear  streamline  coordinate. 

5.  The  boundary  condition  at  the  gas  liquid  interface  was 

suggested  by  G.  1.  Taylor  and  is  Mq  ^ > 

x  3  dfr  ^ ‘f 

since  there  is  no  slipping  (.see  Assumption  .ij . 
u.  The  centrifugal  pressure  gradient  is  small;  thus  =*  O 

dt 


It  should  be  pointed  out  that  the  boundary  layer  equations,  when 
transformed  from  the  x,  y  co-ordinate  system  to  the  f;  and  Yj  system 
(where  £  is  related  to  x  and  ^  is  the  similarity  variable)  have 
neglected  the  partial  derivative  terms  with  respect  to  Similar 
solutions  are  assumed  thereby,  although  strictly  not  allowed,  since 
the  outer  velocity  distribution  over  the  sphere  is  not  of  the 
correct  power- law  form.  It  is  assumed  as  is  usual  that  the  non¬ 
similar  terms  do  not  appreciably  modify  the  results. 


-5 


5.0  FORMULATION  OF  THE  PROBLEM 


Since  the  assumptions  made  correspond  to  reasonable  physical 
circumstances ,  the  mathematical  formulation  for  this  problem 
is  a  straight  forward  application  of  boundary  layer  analysis. 

A  difficulty  that  arises  is  to  determine  what  boundary  conditions 
to  apply  at  the  liquid- liquid  interface  that  describe  the 
physics. 

The  governing  equations  take  the  following  form: 

(A)  Gas  Phase 


(1)  continuity  equation 

1<M>  + 

jT 

(2)  momentum  equation 


(i) 


(2) 


(31  boundary  conditions 


{vo; 


l.  tV  Uc*) 

(B)  Liquid  Phase 

(1)  continuity  equation 


+  MVi  F)  _  o 

(2)  momentum  equation 


(4) 


(5) 
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(3)  boundary  condition* 

At  y  -  0,  the  same  as  gas  phase 

Because  at  the  liauid-liquid  interface,  y  •  -•#  several 
possible  boundary  conditions  can  be  posed;  two  group,  of  them 
are  written  in  following  form. 

Note  that  within  the  "Boundary-layer"  assumptions,  we  assume  that 
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6.0  BOUNDARY  CONDITIONS  AT  THE  LIQUID-LIQUID  INTERFACE 
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shows  the  different  physical  situations  possible  at  liquid- 
liquid  interface.  Consider  Figures  (1)  and  (2)  below. 


Potential  Flow 


Figure  (1)  shows  that  the  velocity  is  zero  at  liquid-liquid 
interface.  Figure  (2)  depicts  the  case  where  the  velocity  is  not 
zero  there.  As  far  as  the  physics  is  concerned  at  liquid- liquid 
interface,  these  are  then  really  two  different  problems.  The 
problem  depicted  in  Figure  (1)  is  a  mass  addition  boundary  layer 
flow  problem,  and  the  one  in  Figure  (2)  is  called  a  "spallation" 
problem.  Therefore,  the  two  different  boundary  conditions  to  be 
specified  at  wj  a  -  °0  are, 

For  problem  (1),  mass  addition:  |  For  problem  (2),  spallation: 


ut=o 

o 

where  K  above  is  a  physical  parameter 
surface  tension  or  viscosity,  both  of 
the  liquid  temperature.  As  the  first 
K  *  constant.  _  g 


0*=  Kw  4^  (6) 
14  =  o  * 

which  can  be  related  to  the 
which  may  be  functions  of 
approximation  we  may  assume 


The  governing  equations  are  just  the  classical  boundary 
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general,  the  standard  non-dimensional  similarity  transformation 
variables  will  be  useful,  since  they  transform  the  two  partial 
differential  equations  into  one  ordinary  differential  equation 
Of  course,  we  must  check  whether  our  boundary  conditions  can  also 
be  transformed  in  similarity  space. 


7.0  TRANSFORMATIONS 


As  is  usual,  we  introduce  the  non-dimensional  stream  func¬ 
tion  ^  ,  so  that  the  gas  phase  dependent  variables  are. 


u  *  A 


vr=  -LdUL 


m 


where  is  the  characteristic  length  of  the  body. 

Defining  the  similarity  variables,  Jp  and  J9  such  that 

r*  i 

%  —  J  ('*  (X)  U  J  y 

1  V-  U- 


(8) 


and 


1 


U(»>  W)A 


(9) 


We  arrive  at  the  standard  stream  function,'  written  as 

r  *  (  U- TJj ) 14  6fo) 

THs  is  therefore  a  definition  of  GC")  . 

All  the  terms  in  momentum  equation  are  transformed  again 


(10) 


in  the  standard  way. 
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U.*  U(y)  G 


(ii) 


IT*  - 


(12) 
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3LgV±  -i 


4it  =  U  r  &  (1) 
^  A(z.- 

Va  StH-L  T1J  k-1  &' 


XJJ  rx 

=  2 .riTT 


’ . 1  Gas  Phase 


Substituting  all  those  _erms  into  gas  phase  momentum  equations, 
an  ordinary  differential  equation  is  obtained  as 

G"  +  G  G*  ♦  if  1)  1'  -  G' l)  =  o  (l#, 

«) 

where 

JS’C* )  as  2.  d 

^  ./  d  ^  (17) 

For  example  at  the  point from  the  leading  edge,  R(zV-  0.235 • 

(See  Appendix  I.) 

The  boundary  conditions  are  now. 

y  =  o  £g-o  *><  q!^\=  ul(x)|  (i8) 

7  —  { c-'w  =  i -0 

Note  that,  in  general,  at  f|  *  0,  Cj>  (fl)  O. 

7.2  Liquid  Phase 

Similarly,  the  non-dimensional  stream  function  and  the  sim¬ 
ilarity  variables  for  liquid  phase  are  defined  as 


(19) 


v-  U5u-7iP  Lty) 

PjiaJiMii. 

J  J»  fit 

7  = 

T Jz^xtf 


(20) 


(21) 


The  only  difference  between  the  above  quantities  equation 
(19-21)  for  the  liquid  and  those  of  the  gas  is  the  kinematic 
viscosity  coefficient  in  equations  (8-10).  Again  a  similar 
procedure  is  used  to  transform  to  the  similarity  space.  The 
governing  differential  equation  is  obtained  for  the  liquid  phase 


as, 


L"  4 


LL*  +#j)  [J[-  -  L'  ]  =  O 


(22) 


Notice  the  term  in  the  above  equation.  The  operator  is 

therefore  different  than  that  of  equation  (16)  for  Gty)  • 
The  boundary  conditions  for  equation  (22)  are  as  follows  For 


(1)  The  mass  addition  problem: 


y-  0  £|_lo 1  =  0;  l'io) -  lf(«)»Gr"(«).r 

=  o 


23) 

(24) 


(2)  The  spallation  problem: 


|:Q  \Jfl)  =  Oj  =(5  3  L  Co)  *  6  M  r  (25) 


12  - 


f2M 


7-**  = 

W 

where  k  Is  the  non-dimensional  form  of  the  parameter  K. 

By  examination  of  both  the  governing  equations  and  boundary 
conditions,  it  is  easy  to  see  that  this  two-phase  boundary- layei 
flow  problem  is  coupled  by  the  interface  conditions.  (See  equa¬ 
tions  (23)  and  (IS)).  Because  our  main  goal  in  this  analysis  is 
to  find  how  much  liquid  has  been  stripped  away  from  the  liquid 
droplet,  the  liquid- liquid  interface  boundary  conditio.:  is  really 
the  most  important  part  in  this  analysis.  If  the  boundary  con¬ 
dition  at  ijs.  —  00  is  defined,  the  subsequent,  interest  is  the 
asymptotic  solution  to  equation  (22)  for  the  liquid  phase,  The 
specification  of  our  imposed  boundary  conditions  there  will  depend 
upon  the  behavior  of  the  asymptotic  solution  of  the  problem. 

In  equation  (23)  or  (25),  we  have  attempted  to  impose  at  l|  *  •  °° 
either  O  or  | L  -  }k  l^^the  asymptotic  solutions  will 
provide  a  check  of  whether  any  or  either  of  these  are  allowed. 


note:  r  ■  p  /v 

g  g 


and  is  typically  less  than  0.1, 
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8.0  THE  PROCEDURE  FOR  FINDING  THE  ASYMPTOTIC  SOLUTION 


One  boundary  condition  we  plan  to  set  at  the  liquid- liquid 
interface  is  that  .  This  implies  that  OO. 

Following  the  usual  procedure  to  find  the  asymptotic  solution 
of  the  governing  equation,  namely 


L'"  *  LL"  ♦  J2(y)[  l'2]  -  0 


we  first  differentiate  the  above  and  obtain 

Lr  i  ll'"  h  L'L"f|  -lj})  =  0 


(27) 


Applying  the  conditions  l'^O  and  to  the 

above  relation  we  see  that 


-+ 


(28) 


Integrating  the  above  once,  we  have 

,  ///  -  CVj 

L  =  e.  '  4 

■ ,u i  > 

Now  since  L.  V“°°)  must  be  zero  due  to  breaking  up  con¬ 
dition  from  the  droplet  at  -  ,  this  implies  that  Bj  »  0  in 

the  above  relation. 


R 


(29) 


Integrating  equation  (29)  once  more,  then 


(30) 
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Also  l"m  must  be  zero  (since  l(-~)=o  );  this  implies 
that  =  0  in  equation  (20'.  The  proceeding  boundary  condition 
means  no  shearing  stress  at^s—  oa,  which  means  that  actually 
some  fluid  has  broken  away  from  the  liquid  droplet 


Continuing,  we  integrate  once  more,  and  find  that 

I 


^  B3 


(31) 


Again  \J (-  ,  the  original  supposition.  Therefore,  =  0. 

Finally  the  solution  for  L(ljj  is 


L  -  ~  t  8^ 

We  note  that  since  C  1  ^en  ®4  3  C'  t^us 

_C.I1 


(32) 


L=  — x.  r  1 

u  c*  <-  -»*  c 

L  is  the  asymptotic  solution  for  the  mass  addition  problem 


(33) 


In  the  process  of  finding  the  asymptotic  solution,  one  can 
make  a  very  important  observation  as  to  the  physical  phenomenon 
occurring  in  the  two  boundary  layer  problem.  That  is,  L"  can  be 
expressed  as, 


L 


// 


(34) 


This  expression  indicates  that  there  is  an  inflection  point 
for  the  velocity  distribution  of  the  liquid  phase  The  following 
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figure  gives  an  indication  of  the  resulting  expected  velocity 
distribution  that  must  result;  if  we  are  to  match  the  boundary 
conditions  at  y  =  0 


The  physics  should  be  carefully  studied  and  investigated, 
and  we  plan  to  do  so.  A  suitable  explanation  for  the  maximum 
velocity  in  the  liquid  layer  is  necessary. 

The  other  boundary  condition  for  the  spallation  problem  we 
plan  to  specify  is  This  implies  that 

l\-«a  *  e ,/A  -f  c  (“ 

The  above  expression  indicates  that  As  ^ 

approaches  ^ approaches  C.  This  corresponds  to  the 

velocity  distribution  in  Figure  2.  Because  this  integ? 
constant  C  is  a  function  of  many  physical  properties,  i. 
difficult  to  determine  it  a-priori.  For  this  reason,  a  f  t 
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study  of  the  asymptotic  solution  of  the, spal lation  problem  is 
required  but  is  not  made  here.  We  plan  to  attempt  such  an 
analysis  later 
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9.0  PROCEDURAL  OUTLINE  FOR  SOLVING  THR  COUPLED  EQUATIONS 


The  equation  of  motion  in  the  gas  and  liquid  phase  have 
been  transformed  into  ordinary  non-linear  differential  equations. 
Because  these  equations  are  notv-linear  and  the  boundary  conditions 
are  given  at  more  than  one  point  (split) ,  the  existence  of  a  solu¬ 
tion  to  this  problem  is  not  guaranteed.  Therefore,  at  this  point, 
an  analytic  approach  which  attempts  to  solve  this  problem  is 
deferred  until  later. 

A  numerical  solution  is,  therefore,  applied  to  this  problem. 
Because  the  boundary  conditions  at  the  gas- liquid  interface  are 
indefinite,  the  procedure  used  in  solving  this  problem  is  to  first 
guess  a  condition  there.  Then  an  iteration  process  is  applied 
until  all  the  other  boundary  conditions  are  satisfied.  For 
example,  at  y  ■  0,  we  guess  the  slipping  velocity,  that  is  a 
value  for  G'(0),  to  meet  the  condition  that  G'  (oo)  ■  1.  A  part 
of  the  solution  is  then  a  value  for  G' ' (0) . 

We  then  proceed  to  the  liquid  phase;  we  start  here  by  setting 
t L{p)*Gf<0)  ,  and  G'(O)  P.  Again  L(»)  -O  .  We 

solve  numerically  the  governing  Falkner-Skan  equation  for  (jjp  , 
and  if  we  are  doing  the  mass  addition  problem,  we  stipulate  that 
aty|m.-<atj  ,  ^(-Oo)csO  “  **  we  cannot  satisfy  this  condition 

at  the  liquid-liquid  interface,  then  we  go  all  the  way  back 
to  the  gas  phase,  guess  a  new  value  for  G* (0) ,  obtain  a  new  value 
for  GM(0),  and  then  solve  the  liquid  phase  with  these  new  con¬ 
ditions.  Eventually  (recall  all  these  are  one  value  of  jkj>)) 
there  will  be  only  one  slipping  velocity  £^(p)  *  ^-(O)  which 
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will  satisfy  both  outer  boundary  conditions,  which  are  Qi'i00)  —  I 
ar.d  .  w*  •'hen  proceed  to  other  points  along  the 

streamline,  where  we  vary  jj  from  the  near  the  stagnation  point 
to  the  maximum  (.shoulder  location)  height  at  */R  -  ir/z 

Mass  Stripped: 


We  can  easily  determine  the  mass  added  to  the  liquid  boundai . 
layer  at  one  J3  location  by  the  relation, 


C363 


Thus  the  total  mass  stripped  in  an  arbitrary  time  period,  T,  is 


_ _  r  06 

^Dj^T  UISSlL  J  dj 


(37; 
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10.0  COMPUTER  PROGRAM  (SEE  APPENDIX  II 


Before  a  numerical  solution  of  the  Falkner-Skan  equation 
can  be  made,  it  must  be  reduced  from  a  split  boundary  value 
problem  to  an  initial  value  problem.  The  boundary  conditions 
for  the  generalized  Falkner-Skan  equation  usually  prescribe 

To  numerically  integrate  the  equation, 
f" (0)  must  be  known.  Note  that  -£Cll)  corresponds  to  either 
or  L(l tj) 

To  solve  for  f"(0),  we  must  use  the  other  boundary  conditions, 
and  a  4th  order  Runge-Kutta  method.  The  following  scheme  is  used. 

where  are  initial  guesses  for  f^'(0)  and  f2"(0)^  j?,  t 

are  the  corresponding  values  of  §  and^^(oo)  The 

^f(o)  so  obtained  is  added  to  $"(0)  .  and  £»{6)  becomes 
yu)  This  process  is  repeated  until  I  f  Ul  -  £1  <  a's 

The  value  of  f"(0)  which  satisfies  this  condition  is  then  the 
correct  value.  This  method,  however,  is  not  full-proof.  This 
entire  process  depends  on  the  initial  guesses  of  £"(0).  This 
initial  guess  must  be  within  50%  of  the  true  value,  or  the 
method  will  not  converge. 


Once  the  method  converges,  a  least  squares  curve  fit  is  used 
to  find  the  relationship  between  f'  '0)  and  f"(0).  The  values 
obtained  are  carved  out  for  a  three  degree  polonomial.  This 
should  be  more  than  sufficient,  since  values  of  f  (0)  are  usually 
less  than  .5.  The  least  squares  curve  fit  was  based  on  a  weight¬ 
ing  function  of  1  for  all  values. 
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Did 

Counter 
<Run  Out 


e  B  30°,  B  ■  . 1850 


A  ■  .67267, 
0 


f'  (0) 

.01 

.02 

03 

.04 

OS 

06 

.07 

.08 

.09 

1.0 

1.2 

1.4 
1.6 
1  8 
2.0 

2.5 
3.0 
3.  S 
4.0 
5.0 


f"  CO) 

.6698 

.6668 

.6638 

.6607 

.6575 

.6542 

.6508 

.6472 

.6436 

,6399 

.6322 

.6241 

.6157 

.6069 

.5978 

.5735 

.5471 

.5188 

.4885 

.4228 


Ax  «  -  0.2777,  A2  ■  -0.50645,  A, 
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»  0.12543 


9  «  6  -  .2950 


f'  (0) 

.01 

.02 

.03 

.04 

.05 

.06 

.07 

.08 

.09 

1.0 

1.2 

1.4 
1.6 
1.8 
2.0 
2.  5 
3.0 

3.5 
4.0 
5.0 

Ao  ■  .77061,  A1  *  -  0.38473, 


f'(0) 

.7666 
.7627 
.  7586 
.  7544 
.7502 
.7458 
.7413 
.7368 
.7327 
.7274 
.7177 
.7076 
.6972 
.6864 
.6753 
.6460 
.6148 
.5816 
.5466 
.4713 

A2  ■  -  0.47954,  A3  »  0.10396 
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R  ■  0.0808 


A  =*  .56626 


a  .  e.n°  t 


f  CO) 


f"(0) 


01 

.5646 

.02 

.5630 

.03 

.5613 

.04 

.5594 

.05 

.5575 

.06 

.5554 

.07 

.5522 

.08 

.5509 

.09 

.5485 

1.0 

.5460 

1.2 

,5407 

1.4 

.5349 

1.6 

.5288 

1.8 

.5223 

2.0 

.5154 

2.5 

.4966 

3.0 

.4756 

3.  5 

.4526 

4  0 

.4276 

5.0 

.3721 

A.  *  -  0 . 14753,  A2  “  ‘  °-565176»  A 


«  . 167826 
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S  ■  90°.  8  ■  0.0 


A  -  .46984, 
0 


f'  (0) 

.01 

.02 

.03 

.04 

.05 

.06 

.07 

.08 

.09 

1.0 

1.2 

1.4 
1.6 
1.8 
2.0 

2.5 
3.0 

3.5 
4.0 
5.0 

A1  -  -  0.0095135, 


f"(0) 

.4695 

.4603 

.4689 

.4684 

.4677 

.4669 

.  4660 

.4650 

.4638 

.4625 

.4595 

.4561 

.4522 

.4479 

.4431 

.4295 

.4135 

.3953 

.3750 

.3287 

■  -  0.66554,  A 


-  .24078 
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This  report  has  dealt  with  the  problem  of  the  interaction 
between  a  liquid  droplet  and  the  convective  flow  field  surrounding 
it  This  problem  has  received  some  attention  by  other  investi¬ 
gators,  as  indicated  in  the  list  of  references.  Most  studies 
have  been  experimental  ones,  and  one  in  particular,  by  Nicholls 
and  Ranger,  has  attempted  an  approximate  boundary- layer  analysis 
by  first  assuming  arbitrary  simple  velocity  distributions,  both 
in  the  gas  and  liquid  layers. 

This  report  is  an  attempt  to  solve  the  problem  of  aerodynamic 
stripping  of  liquid  from  a  spherical  drop  by  solving  (numerically) 
the  full  coupled  steady  boundary  layer  equations  (for  incompressible 
flow)  and  stipulating  certain  realistic  interface  conditions, 

No  solutions  are  given  since  the  work  represents  a  funded  effort 
of  only  one  half  year.  However,  considerable  insight  into  the 
nature  of  the  coupling  behavior  of  the  gas  and  liquid  viscous 
flows  has  been  made,  and  we  expect  later  to  carry  out  to  com¬ 
pletion  the  solution  to  the  model  proposed.  The  main  parameter 
of  interest  to  be  solved  for  will  be  mass  of  fluid  in  the  cir¬ 
cumferential  liquid  layer  that  was  swept  along  by  the  gas  stream. 
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APPENDIX  I 

CALCULATION  OF  THE  VALUES  FOR  S 


X,,  -  2jU.  £  M.  -L 


V,  -  >  U’ 

V-  =  R  S\H  f 1 


q  ^rZ  ..  J.  X 

TIL  *i«  ^ 
R1* 


dXJi 


ss  2.  u<* 


cos 


t  ,  f*  Rz»»*  -fc  •  i s"l~p~ 

^  ^  ru« 

Jo 

-  f-p-  “■♦•‘HO 


■$c 


f’  +  fill 


-  COS 


il 
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2T 


J  a?  4 *• 
2-JU ,  7r(i:^c°s'R^  fC§ 

-  4i.  ^  Co*  ^  _A - 

-  3  J  *„».$.  lrl  1^ 


At 


TL 

2- 


£l£l  *  0 


j: 

TL.  0.o8\ 


oouod  !  j  dooooooooooo 


COMPUTER  PROGRAM  FOR  SOLVING  THE  FALKNER-SKAN  FOUATIUN 


_ IMPLICIT  RFAI  gSIA-H.O-YI _ ; _ 

DIMENSION  Y1 ( 200 ) ,Y2<  200 ) , Y3 ( 200 ) ,DY 1 ( 2  00 ) , 0Y2 I ?00 ) ,0Y3{ 200) , XING  I 
12  00  )  .ALPHA  I  10)  iflETAAUO  ) ,  TR  Y  ( IQ  I  .  THETA2  I  1 00  ) ,  ThF  T  A1  ( 1  00  )  ,  PR  A  I  20  ) «  L 
1TX( 2) .ZTYI2) ,2F1 ( 200) ,Z F? ( 200 > . ZF3 < 200 ) ,ZFa< 200) ,ZXIMC(20fl ) 

CALL  UNDERZ ( 'OFF* ) 

CALL  ERRSFT ( 207.256.-1 ,1 ) 

L _ _ _  ..  .  ... 

t 

INHIALJZF  VARIABLES 


READ! 5,10) ( PR A ( J ) , J=1 ,20 ) 

10  FORMAT (20F4. 2) 

H=5.0D-2 
DNSTY=1.D-1 
DNST  Y= 1 .000 

_ TRY  (  1  1=0.000 _ _ _ _ 

TR Y ( 2 ) =0 . 29500 
TRY ( 3 ) =0.080800 
TRY<4)=0. 10500 
ZF2  (1)*9Q. 

Z  F2 ( 2 ) =45  . 

_ _JLE2jLiL=iia, _ _  _ _ _ 

Z  F2 ( 4 ) =30 . 

X  =  0  «.QD.O  .  .  .  ... 

M  =  1  01 

_BETA=1.QQQ  . 

no  5000  J K L  =  1  ,4 

_ Bl£IA  =  TRY  l.JKl_) - - 

WR I TF ( 6 , 5 06 )  ZF2 ( JKL ) , TR Y ( JKL  ) 

506  FORMAT (' 1 ' ,23X, 1  THETA  =  •  .  F 5 . 1 .  •  DEGREES., — BETA  =  ■  ,  F 7 . 4 ,  /  ,  /  ) 
WRITE (6,507) 

507  FQRMAI.U _ '.,23X,  '  F  '  '  (  0)  «  ,  10X,  '  F  "  ".(  0 )  '  ,  /  ) 

DO  4999  LM= 1,20 


NOTF:  THE  VALUE  OF  CONST  IS  TO  BE  IJSFD  AS  FOLLOWS  :**##*##**#>»w 

IF  BETA=0.0  THEN  C0NST=1.0  OR  0.5  DEPENDING  ON  THE  FORM 

_ OF  THE.. BIAS! US  ED  1.1  A T  I  ON  TO  BF  lISFn - - - - - - 

IF  BETA  0.0  THEN  C0NST=1.0  ALWAYS!!!!! 


C  ONS  Ta  5. 00^  1 

CON  ST  =  1 . ODO 
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FINFKNaO.ODO 

- s-mfeiuui.oao - 

Y 1  (  1  )»0.000 
Y2iiua.aino 

Y  2 (  1  )“PRA( LM) 

ZFULM>»PRA1LM>  - 

AL  PHA ( 1 J-0.5D0 

- AU»-WA.i2  )jtft^fcDO _ 

INCRFM*0 

— icotwua — -  - - 

c 

c 

c 

— C - maim _ program  *************** 

C 

c  .  _ _  . 

c 

HOLD 1 “ALPHA  1 1 )  _ _ _ 

HOLD2“ALPHA( 2 ) 

_ A52A _ QJD  hid  LLal  >26.2 _ 

IF (INCREM.EQ.2)  GG  TO  6314 

. .  .  .GO  ...TH.  A3 15 _ _ _  _ 

6314  ALPHA( 1 )“H0LDl+3.0D-2 

AlP.MA(2)  =  HOLD2  +  3.QD-2  . .  .  _ 

H0LD1“ALPHA< 1 > 

_ HHI  H?rAI  PHAi  ?  > _ 

INCREMaO 

. . .  . 1  COUNT.*  Q. _ _ _ _ _ _ _ 

6315  DO  150  J J  =  1 1  2 

Y3U  1*ALPHA(JJ) _ _ _  ..  .  . 

KKK  = J  J 

_ Xi.NC.llJ.~X _ _ _ _ 

DY1  ( 1 )«Y2<  U 

_ nyz(  i  )cV3<  1 ) _ _ _ _ _ _ _ _ 

DY3< l )«-C0NST*Yl(l)*Y3(  1 )- (  BETA*  (  DNiiT Y- (  Y2  (  l  )**2)  )  ) 

_ ,DQ_X_Jf.2-*M _ - _ _ - - -  - 

KN»J 

- PlsHADY.lCJ-1  ) - 

P2  =  H*DY2<  J-l ) 

_ B3.aH*ay.3  lJ.-r.li. -  -  . 

A*Y1 (J— 1 )+ ( 5.0D-1 )*P1 

_ B-slYZI  J- 1  1+  L5.*00=a  l-*£2. _ _ _ _ _ 

C=Y3(J-1J+(5.0D-1)*P3 

_ QlsHtB _ 

Q2=H*C 

_ _ Q3=H*XL-CQNSI1*A*C-J  B ETA*  (  ONSTY-I  B*.*2J-UJ . . 

A=Y1 ( J-l >+<5.0D-l >*Q1 

B  =  Y2  l  J-l  1  ♦  1  5 »jQ.D."l 1*02 _  _ _ _ 

C=Y3(J-1)+(5.0D-1)*03 

_ aisHSB _ 

R2=H*C 

_ RJjgH^i  l  -CQNS  r.l^AYC-l  BETA*  1  DMSXY-.LB * *2 1.1  U _  _ 

A  =  Y1 (J-l ) +R 1 

_ BaY2XJ.-rl.li.R2 _ _ _ _  - 


C  =  Y3( J-l ) +R  3 

S1*H*R _ 

S2=H*C 


32 


ISO 


MO 

C 


1 


SS»H*( <-CONST>*A*C-(BFTA#<DNSTY-<B**2) ) ) ) 
Y  1  i.iicyi  li.  l  ujoun  ODOJJi’Gi- 


1 

3000 


.3ao- 


Y2< J)=Y21 J-1)+!P2+ (2.000) *02+ (2.000) *R2+S2 1/6.000 
Y3(  J)  =  Y3(  J-l)+(P3+(2.QDO)*Q3+(2,QD0)*R3+S3)/6.ODO 
n V1 ( J  ) =  Y2 1 J  I 
nv2(J)=Y3(j) 

DY3(  J»»-C0NST*Y1  (  J)+Y3(  J  ) - ( BFTA* ( DNSTY- I  Y2  (  J)**2>  )  ) 

.,1L1Y21.JU£T,.1.5DQ)  SO-TP  .3000 _ 

XINC! J)=XIMCI J-l )+H 
CONTINUE 

IFIKKK.RO. 2)  GO  TO  390 

FINF1=Y2(KN) 

on  to  iso 

.P-LNE2_TY2i_*NJ _ _ 


C0NT  I  NIIF 
XB=FINFKN-FINF2 

XC=(  ALPHA!  2) -ALPHa  (  1  |  )  /  (  F  I NF  2.-F  I  NF  1  ) 
PELTF=XC*XB 
ALPHA! 1 1  =  AL  PHA ( 2 ) 

ALPHA(2)=ALPHA(2)+DELTF _ 

1 F( OABS ( F INFKN-F INF2 I .LE .1 .00-5 )  00  TO  SOS 

I  COUNT  = 1C0UNT+1 . 

IF!  I  COUNT. EO, 10)  I  NCR  FM  =  2 
CONTINUE 


PRINT  OUT  RESULTS 


JillS _ WR I T  F ! 6 . 5 04  )  Y2(l).Y3ll) 


504 

4999 


SOOO 


FORMAT! '  S23XfF14.10t5X.F14. 10) 

.CONTINUE  ....  .  . . . 

I  COUNT  =  0 

-JNCREHsQ . . 

CONTINUE 

-S.IXLE _ 


END 

/.*  .  . . 

//go.sysin  no  * 

.01  .02  .03  . 0.4 _ .  05  .  .06  .07  .08  .09  .10  .12.  .14  .16  .18  .20  .25  .  30  .35  .  4. 
/* 
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